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Abstract 

The CP violating parameter e'/e is computed using the low-energy dynamics of the chiral theory 
supplemented by vector resonances. The divergent contributions coming from strong ir-ir scattering 
are tamed by vector- meson exchange terms. This amounts to softening the fast growing high-energy 
behaviour of tt-tt scattering. The final result for e'/e shows a smooth dependence on the cut-off 
where low energy dynamics is matched with that of QCD. 

PACS numbers: 13.25.Es, ll.15.Pg, 12.39.Fe 

Keywords: e'/e, Direct CP- Violation, K meson, Large- N c QCD 



* Email akundu@cucc.emet. in 

t Email: paschos@physik . uni-dort mund . dc 

■f Email: zhou@zylon . uni-dort mund . de 



1 



1. Introduction. The decays K — > tytt are best described by a low energy effective 
Hamiltonian 

H = ^ {|>(?V) + rw(?V))Qi} (1) 

with Zi(q 2 ,fi 2 ) and 2/i(<? 2 ,/i 2 ) being the Wilson coefficients and £ g = V* s V q d, t = £t/£« ■ Qls 
are 4-quark operators. For the definition of the operators and other notations, see ref. ffl 
which we closely follow. Matrix elements for two of these operators, Q 6 and Q 8 , are most 
important for the evaluation of e'/e: 

Qe = -2 s(l + 75)g?(l-7sK 

q=u,d,s 

Q* = - 3 e q s(l + l5 )qq(l- lB )d, (2) 

q=u,d,s 

where e q = (2/3,-1/3,-1/3). The QCD corrections included in the Wilson coeffi- 
cients represent the short distance terms computed in perturbative QCD. They depend 
on [ln(Q 2 //i 2 )] 7///3 and to next-to-leading order (NLO) corrections in a more complicated 
way. The numerical values have been tabulated by various groups 0, HI . Comparisons of 
the results show that the various groups agree with each other but values for the coefficients 
depend on the renormalization scheme. The //-dependence in the coefficients is expected 
to be cancelled by the scale dependence of the matrix elements of the operators introduced 
through the upper cut-off in the integrals, and the running strange quark mass. 

The matrix elements of the form (7T7t|Q6|A') an d (7T7T | Qg | i^) include tree level contribu- 
tions and loop corrections. These are low energy processes which must be dealt with by 
methods other than QCD. Our method is to use the low energy chiral theory for calculating 
tree and loop diagrams and then match the results with the short distance contribution, i.e. 
the QCD scale /i is matched with the upper cut-off A c appearing in the chiral loops. An 
important criterion for the success of the calculation is smooth (and weak) dependence of 
the results on \i = A c . 

In the large N c approach factorizable and non-factorizable amplitudes are treated sepa- 
rately 0] with the factorizable amplitudes defining the renormalized coupling constants. In 
a Dortmund-Fermilab collaboration 0], it was shown that to O(p /N c ) the divergences in 
the matrix elements of the Q$ and Q$ operators are logarithmic and occur in nonfactorizable 
diagrams. 

The numerical results of this approach at O(p /N c ) were presented in table I of ref.pj, 
which we also adopt in the present article. The results of the diagrammatic method were 
reproduced in the background-field method 0. Let us denote the nonet of pseudoscalar 
meson by the matrix II = P a A a , where A a 's are the usual Gell-Mann matrices; then it was 
shown that to 0(p°/N c ) 

/ " F n an / " F K U 

where (ttq, K ) and / are the bare pion and kaon fields and decay constants, while (n r , K r ) 
and F w , F K are renormalized fields and decay constants, respectively. 
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A large correction in the earlier calculation [5( originates from rescattering of the pions, 
i.e., K — > mi — > mi 1 where the first step involves the weak operators Qq or Q$ to O(p 2 /N ) 




FIG. 1: Feynman diagram for K — ► 7T7t with strong final state interactions. 

and the second process is the strong pion-pion scattering as shown in Fig. ^ The large 
dependence of the cut-off resides on the contact ti — n scattering which is known to have 
a bad high-energy behaviour violating unitarity and needs to be moderated by some other 
amplitudes which restore unitarity. 




FIG. 2: Feynman diagram for K — > tttt with a vector-meson exchange. 

A standard prescription to restore unitarity is to introduce vector-meson exchange dia- 
grams. For the titi — > titi scattering we shall use the contact and the p exchange diagrams. 



The initial state interactions are expected to give smaller contributions, which we will present in a future 
publication Q 
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We accomplish this by using a chiral Lagrangian for pseudo scalars and enlarged by the 
introduction of vector mesons 0, 0, . We extend the calculation of the one-loop diagrams 
with a strong vertex with the addition of a p-exchange diagram. The p is included to repre- 
sent the effects of even heavier vector mesons (like K*). In addition the pions are in I = 
or / = 2 states and the exchange of p-mesons appears only in the t-channel, see FigEJ 

In order to restore unitarity we shall demand that quadratic divergences cancel between 
the contact and the p-exchange diagrams. It is indeed heartening to note that they come 
with opposite signs, and cancel exactly if the following relation is satisfied 



Here h is the pirn coupling strength and / is the pion decay constant (~ 92MeV). The 
logarithmic divergences still remain and should be matched to the QCD logarithms. This is 
our proposal for moderating the high energy growth of 7r — tt scattering. 

Thus, we calculate the one-loop K —> nn amplitudes with both contact and p-exchange 
diagrams, demanding that the quadratic divergences cancel betwen these two sets. The 
value of h ~ 4.8 obtained from Eq.(0J) is slightly smaller than the one obtained from the p 
decay width, but remember that p is only a symbolic representation of all possible vector 
resonances. Since only logarithmic divergences will be present in the final result, the varia- 
tion of e'/e with the cut-off A is expected to be weak. As the weak vertex (with Q Q or Q 8 ) 
is common to both the contact and the p-exchange diagrams, the cancellation of quadratic 
divergences is respected by both operators. 

2. Framework. The effective Lagrangian for pseudoscalar mesons relevant for K — > titt 
decay up to 0{p 4 ) is given by [Io| : 

C eff = £ (id^d^U) + |^(ln[/ f - hif/) 2 + r(MU* + UM)^j + r 2 H 2 (M 2 ) 

+rL 5 (d lx tfd fl U (MU + U ] M)) + r 2 L 8 (MUMU + MU ] MU ] ), (5) 

with (A) denoting the trace of A and M = diag(m u , m<2, m s ), f and r are free parameters 
related to the pion decay constant F n and to the quark condensate respectively, with r = 
-2(qq)/p. 

The matrix U is given by 

U = exp{m/f), (6) 
where the pseudoscalar meson nonet II is given by 

At + -La V + &brf V2n + ^2K+ \ 

n = A a P a = V2n- V -tt° + j^arj + ^brf V2K° (7) 

where A's are the usual Gell-Mann matrices, P a are the pseudoscalar fields, and 

a = cos 9 — V^sinO, b = —= sin 9 + cos#, (8) 

V2 
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9 being the r\ — rf mixing angle. We include i] and 77' contributions in our calculation. It 
is easy to see that though the operator Qq vanishes at tree-level due to the unitarity of U, 
it still has nonzero contributions at the O(p /Nc) level. The loop expansion of the matrix 
elements is a series in 1/ f 2 ~ 1/N C , which follows from the short- distance expansion in terms 
of « s /vr ~ 1/N C . 

There have been numerous calculations of e'/e which try to improve various steps 
The expression of e'/e can be written in a compact notation as 



e' G F oj 

-Im& 



n - -n 2 



[u = 1/22) (9) 



e 2 |e|ReA) 
with 

Ho = [ y^l/i(^)(Qi)ol (1 - toy+rf) » 

i 

n 2 = \ J2yi(v)(Qihl (10) 

i 

The isospin breaking effect (m u 7^ rrid) is taken into account by £I v +t/- 

Our aim is to introduce vector mesons in terms of a Lagrangian which satisfies the low 
energy current algebra. One consistent method is in terms of a non-linear chiral Lagrangian 
with a hidden local symmetry j|6]. In this theory the vector mesons emerge as dynamical 
vector mesons. The three point vector-pseudo scalar interaction is given by 

jiV^Pd^P-d^PP)), (11) 

where h stands for the vector-pseudoscalar coupling. Some typical vertices of p's to pseu- 
doscalar mesons are 

n + (p 1 )n-(p 2 )p° : h( Pl -p 2 )^ 
7i + (p 1 )7i (p 2 )p- : h{p l -p 2 )^ 

K + ( Pl )K°(p 2 )p- : A(p x -p 2 )^, etc, (12) 

which is directly related to the p decay width: T(p) = /i 2 (|p 7r |) 3 /(67rmp), where p,,- is the 
momentum of final state pions in the p rest frame. With T(p) = 149.2 MeV, we find h = 5.95. 
We note in passing that the Kawarabayashi-Suzuki-Riazuddin-Fayyazuddin relation gives 
the value h = ra p j (-v/2L)(l3|. Thus the value of h in Eq.(jU and the two values in this 
paragraph differ by small amounts (~ 19%). The strong four-point vertices involving pions 
are obtained from the first two terms of Eq. (|5|i. The weak vertices are obtained from the 
definitions of Qq and Q$. In the numerical work we shall use the value of h from Eq.(J3J) and 
also h = 5.95 obtained from the decay width. 

We repeated the renormalization procedure and found the following results. For the self 
energies of the pseudoscalars, momentum independent terms combine with the bare masses 
to define the physical masses. A momentum dependent term is included in the wave function 
renormalization and is the same for it and K. The renormalization of F n and Fk is the same 
as in ref. [l|, i.e. there is no h 2 contribution, which leads to the same value for L5, similarly 
the value for L5 — 2L§ is again very small. The quadratic divergences of the factorizable 
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diagrams for (Qe)o, {Q§)2 and (7r°7r°|Q8|-^) cancel out, what remains of them are small 
corrections because to O(p ) these matrix elements vanish. The quadratic divergence from 
the factorizable diagrams of (tt + it~\Qs\K) cancel against the corresponding diagrams with 
vector meson exchanges when we invoke the condition in Eq.(j3J). The surviving term is small 
in comparison with the O(p ) contribution of (7r + 7r — |Qs 1^) - 
We use the following numerical inputs: 

m n = 0.137GeV, m K = 0.495GeV, m p = 0.771GeV, 

7 = ^ = 0.0924, m s (m c ) = 0.115GeV, a s {m z ) = 0.117. (13) 

The strange quark mass has an error of 0.020 GeV [13]. The average quark mass rh is given by 
m = m s /24.4. We also use L B = 2.07X10" 3 , L 8 = 1.09xl0~ 3 , Jm(&) = (1.31±0.10) x 10" 4 Q 
and the isospin breaking factor of £l v + v ' = 0.15 

One can extract ^q CD from as{m z ) at either the continuum upper limit [l5| (m b (rfi b ) = 
4.5 GeV, m c (m c ) = 1.4 GeV) or the continuum lower limit (m^m^) = 4.0 GeV, m c (m c ) = 1.0 
GeV): 

A^ CD = 0.279 ± 0.029(upper limit), ^qcd = °- 275 ± 0.029(lower limit). (14) 

We take, as a conservative estimate, J^q C d = 0.277 ± 0.031 GeV {i.e., between 0.246 and 
0.308 GeV). 

The Wilson coefficients were tabulated |5( for various renormalization schemes and the 
values of A.q CD as functions of the renormalization scale /x. The values show a convergence 
among the schemes as \x increases and approaches the value of \i = 1 GeV. This is as expected 
since QCD is valid at higher momenta. 

A second issue is the matching of the coefficients in the various schemes to the cut-off 
scale of chiral theory. A method for relating the two scales was suggested in The 
method introduces 

1 = ~2 2 2 2 ( 15 J 

and uses the first term as the infrared regulator of QCD and the second term as the cut- 
off for the chiral theory. This approach provides a matching of the two scales A c and fi. 
Recalculation of the evolution of the coefficients 0| brings the values of the HV scheme 
closer to NDR, which are anyway close to the leading order results. All this motivates us to 
use the values of the NDR scheme. We shall use values for Aq^ d = 0.245 GeV, however, we 

check that interpolation to J^q CD = 0.277 ± 0.031 GeV changes the values of e'/e at most 
8%. Althernative ways for matching the two theories have also been introduced in other 
articles [l7| . 

3. Results. As mentioned already, a previous work demonstrated that renormalization 
of physical quantities (wave functions, masses and decay constants) render the factorizable 
contribution to (Qe)o,2 and (Q 8)0,2 to O(p /Nc) finite. There are loop corrections introduced 
by the non-factorizable diagrams which to order p°/Nc were found to be logarithmic. Going 
one step further corrections of order p 2 /N c were studied 0, arising from the contact terms 
which have a quadratic dependence on the cut-off scale A^. We combine the contact terms 
with the vector meson exchange diagrams and cancel the quadratic divergence. 
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We present in this section the results for the contact terms and vector meson exchange 
diagrams to order p 2 /Nc in terms of integrals which are summarized in Appendix^! in order 
to make the reading easier we give in the text explicit formulas for the decay K° — > 7r°7r° 
where the results are shorter. For the decay of K° — > 7r + 7r~ we collected the results in 
Appendix [BJ In both reactions we included the n + ir~ and 7r°7r° intermediate states. 

The contact terms for K°(pk) — > 7r + 7r~ — > 7r°7r° give 

2r 2 

iM^x = i^-^[AI 9 (m n ,m^,p K ,p K ) + BI u (m n ,m n ,p K ,p K ,p K ) 

+ACIs(m n ,m n ,p K ) + BCIg(m n ,m n ,p K ,p K )} (16) 

with A = -%L 5 m 2 K , B = 8L 5 , C = (xi + X2)/4 + m\ - m 2 and Xi = rm i- 
The contact term for K°{px) — > vr 7r° — > 7r°7r° is 

r 2 

^con2 = i ^^j. 3 C'[AI 8 (m n ,m n ,p K ) + BIg(m w ,m w ,p K ,p K )] (17) 

with C = (X1+X2) • The functions Ii(rrij, rrik,p, ■ ■ ■ ) etc represent four dimensional integrals 
which we define in the Appendix El The notation with the numbers as subscripts follow 
the convention introduced in two Ph.D. theses at Dortmund University 18], where explicit 
formulas for the functional forms after integration are included. 
The p— exchange diagram for K°(pk) — >• 7r + 7r~ — > 7r°7r° is 

2h 2 r 2 f 1 
+AI s (m w ,m p ,p 1 ) + BI 9 (m 7r ,m p ,p 1 ,p K ) 

+2AI 30 (m n , m n , m p , p K , p u p x ) + 2.6/3! m n , m p , p K , p u p K , p x ) 
+2(m 2 K - m 2 r )[AJ 2 9(m 7r ,m 7r ,m p ,p^,pi) + 5J 30 (m 7r , m n , m p ,p K ,Pi,PK)]} 

(18) 

Finally the p— exchange diagram for K°(pk) — > tt tt — > 7r 7r is zero 

^exch 2 = (19) 

because the 7r°7r°p vertex does not exist. 

Including the vector mesons with the condition in Eq.(j3J) eliminates the quadratic depen- 
dence on the cut-off. This is our method for regularizing the integrals in terms of physical 
particles and interactions which preserve the symmetries. The remaining logarithmic de- 
pendence of the cut-off will be matched with the ln/z dependence of the QCD. 

We give in table |H the contributions to 0(p 2 /Nc) from the contact and the p exchange 
terms for (Qq)o and (Qs)2 i n un h of r 2 ■ MeV as a function of A^ in the interval Ac = 0.7 
GeV to Ac =1.0 GeV. The cut-off scale must be larger than the mass of p and the first 
column is given only as a point of reference. We note that the dependence of (Qe)o um and 
(Qd)T m 011 Ac is very small. Since the value of h from Eq.(JJJ) is smaller than the value 
obtained from the p — ► tttt decay width, we repeated the calculation for h = 5.95 in table 
ITTl corresponding to the coupling from p decays. The values for e'/e are slightly smaller and 
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A c — U. / (je V 


A c — u.o (je\ 


A no r^aVA 

a c — u.y (je\ 


a in r^oAY 

A c — l.u (je\ 


• i n \ con 
1 \V6/ 


— 14. o 


17 C 

— 1 1.0 


on a 


— Zo.4 


„' /n \p 
1 \V6/ 


O.D 


o.y 


11.0 


14.0 


4 in \ sum 


Q Q 
— O.O 


— o.O 


Q Q 
— O.O 


C Q 
— O.O 












i /n„\ con 

* \V8/2 


0.Z4 


1 .4o 


a 7 

O. ( 


i n i 

1U. 1 


i(Q*) p 2 


-2.30 


-3.15 


-4.11 


-5.17 


i(Qs)f m 


3.94 


4.28 


4.59 


4.93 












Total e'/e(10- 3 ) 


2.23 


1.84 


1.53 


1.2 



TABLE I: The contact term and the p— exchange contributions to 0(p 2 /N c ) for the matrix elements 
of (Qq) and (Qs) (in units of r 2 • MeV) as well as e'/e as functions of the cut-off scales A c . The 
value of h is taken from the cancellation condition of Eq.@ 





A c = 0.7 GeV 


A c = 0.8 GeV 


A c = 0.9 GeV 


A c = 1.0 GeV 


• //O \ con 


-14.8 


-17.5 


-20.4 


-23.4 


i(Qo) P o 


9.93 


13.6 


17.7 


22.3 


■ 1 r\ \ sum 


-5.36 


-3.9 


-2.7 


-1.1 












■ //n \ con 


6.24 


7.43 


8.7 


10.1 




-3.51 


-4.80 


-6.27 


-7.89 


i (Qs) s 2 um 


2.73 


2.63 


2.43 


2.21 












Total e7e(l(r 3 ) 


2.03 


1.57 


1.19 


0.8 



TABLE II: The contact term and the p— exchange contributions to 0{p 2 /N c ) for the matrix ele- 
ments of (Qe) and (Qs) (in units of r 2 • MeV) as well as e'/e as functions of the cut-off scales A c . 
The value of h is taken to be the physical one h = 5.95. 



the variation of the matrix elements with the cut-off is larger. For the calculation of e'/e 
we use, for the tree and factorizable contributions the values from table I of ref. 0, which 
are primarily responsible for the remaining A c dependence of e'/e. The results reported 
in this article present a complete calculation of the matrix elements Qq and Qs to order 
p 2 /Nc- The presence of the vector mesons restores to a large extent the unitarity of the 
theory and acts as an upper cut-off for the integrals. Our results suggest that a non-linear 
chiral lagrangian with a hidden local symmetry may be a more suitable low energy limit for 
QCD. 

As mentioned already, the values of the matrix elements are very stable. The calcula- 
tion of e'/e uses the coefficient functions of NDR at A.q CD = 0.245GeV and m s (lGeV) = 
0.125GeV. We found an imp roved stability of the values for e'/e which are consistent with 
the experimental results JjJ, |20( • The main conclusion is that the presence of vector mesons 
improves the calculation of the matrix elements by making them more stable functions of 
the cut-off. 
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We demonstrated that the chiral theory enlarged by the introduction of vector mesons can 
eliminate quadratic divergences to 0(p 2 /N c ). The improved stability of e'/e is encouraging 
to extent the calculation to the initial state interactions. We expect the changes to be small, 
but we plan to complete them and present them in a longer article [9]. The extension of 
the method to the amplitudes A and A 2 will involve additional operators Q\, Q2, ■ ■ ■ with 
considerable increase in the computational work. It will be interesting, however, to find out 
whether vector mesons make these amplitudes also more stable. 
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APPENDIX A: FOUR DIMENSIONAL INTEGRALS 

Several integrals have been used in this article and we try to define then in a compact 
notation. The integrals I3, J4 have the same denominator but have different numerators 
separated from each other with semicolons 



,4 {i;(p-g)> 

' 3;4 ~ (2vr) 4 J aq {q-ky-m* 



d % M2 2 ( Al ) 



The integrals Is, Ig, ho, hi and lyi have again the same denominator but have different 
numerators separated from each other with semicolons 

r i f ,4 {lj (p ■ g); ff 2 ; (gi ■ gXgg • g);g 2 (p • g)} ( , , 

-<8;9;10;11;12 — (n U \ « 9 7 9 oTF? Tv? 21 \ AZ ) 



(2tt) 4 J H (g 2 - ml) [(q - k) 2 - m 2 2 ] 
The same notation is used in the integrals I29, ho and /31, 

j i f ■a {l;(Pi-g);(Pi-g)(P2-g)> / a o\ 

^29;30;31 J « Q ^ _ ^ [(g _ fc )2 _ m |] [(g _ p )2 _ m |] 
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Among these integrals I3, I4, hojn and I± 2 have quadratic divergences in the cut-off regu- 
larization scheme. The quadraticlly divergent parts are given by 

I 3 (m, k)\ A 2 div = y^^Al 



(4tt; 



h(m,k,p)\ A 2 div = ^^ A 2 



2(4vr; 



ho{m 1 ,m 2 ,k)\ K 2 div = j^-^ A 2 



(4vr; 



I 11 (m 1 ,m 2 ,k,p 1 ,p 2 )\ A 2 div = T^j^ Ag 



4(4tt) 
(k ■ p) 

I 12 (m 1 ,m 2 ,k,p)\ A 2 div = ——A 2 C (A4) 



2(4tt) 



Using the quadratic divergences and the formulas in this article the reader can verify the 
cancellations. 



APPENDIX B: K -> Tr+vr" DECAY AMPLITUDES AT 0{p 2 /N c ) 
The contact term for K°(p K ) — > 7r + 7r~ — > n + (pi)n~ (p 2 ) is given by 



^2 



j^li = ~^ 3v /- 3 [^10 ("t»r, rn 7V ,p K ) + BI 12 (m 7V ,m 7r ,p K ,p K ) 

- 2AI 9 (m n ,m n ,p K ,2p 2 - p^ - 2BI u (m n , m n ,p K ,p K 2p 2 - p x ) 

- ACI s (m n ,m n ,p K ) - BCI 9 (m n ,m n ,p K ,p K )] (Bl) 



with C =(xi + X2) + (mi 

The contact term for K°(p K ) — > 7r°7r° — > 7r + (pi)7r~(p 2 ) is 
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«M c + on2 = ~» ~ [^10 (^vr , m w , Pif ) + BI^TTl*, 777^, p^, p K ) 

- ,4/9(777^, m^, Pa:) - BI u (m n ,m n ,p K ,p Kj p K ) 

- AC'I 8 (m n ,m n ,p K ) - BC'Ig(m n ,m n ,p K ,PK)] (B2) 

The p— exchange diagram through 7r + 7r~ loop gives 
/i 2 r 2 f 1 

^c + xchl = (~^^7fJ |-^2-[^3(^p,Pl) +5/ 4 (777 p ,p 1 ,p^)] 
+ Al^m^m^p!) + 73/9(777^, 777,,, pi,p*;) 

+ 2AI 30 (m n ,m 7r ,m p ,p K ,p 1 ,pi) + 2BI 31 (m n ,m 7r ,m p ,p K ,p 1 ,p K ,Pi) 

+2(m 2 K - ml) [A/29 ("7tt, "V, m p ,p^,pi) + BI 30 (m n , m n , m p ,p K ,Pi,PK)]} (B3) 

with C" = (xi + X 2 )/4 + (777^ - m 2 ). 

The p— exchange diagram through 7r°7r° loop gives the same contribution, i.e. 

-W c + xdl2 = AC; hl (B4) 
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It is straight forward to verify that the cancellation condition of Eq.(j3J also holds for 

K — > 7T + 7T~ . 
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